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Abstract 

We investigate the NBI matrix model with the potential XA + X^ 1 + (2n + 1) log X 
recently proposed to describe IIB superstrings. With the proper normalization, using 
Virasoro constraints, we prove the equivalence of this model and the Kontsevich matrix 
model for n ^ and find the explicit transformation between the two models. 
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1 Introduction 



The investigation of matrix models with an external field and a logarithmic potentials was 
initiated in [|T], 0. These matrix models are related to the so-called NBI matrix model, which 
appeared recently in the context of the IIB superstring matrix model |3|, [|, |5|, [J. In |5|, 
the following (M)atrix model action was proposed: 

%BI = ~\ tr Y- 1 ^, A v f + (3 tr Y + 2rj tr logF - X - tr ttP 4 ^, *]. (1.1) 

It possesses M = 2 supersymmetry in the large N limit || . The matrix Y plays here the role 
of the world-sheet metric to be integrated out in order to obtain the effective action. Then, 
the (nonlocal) logarithmic term is a curvature term of the world-sheet metric. As was shown 
in M , the relevant choice of the constant r] in front of this term that leads to a non-Abelian 
Born-Infeld action for the string coordinates is rj = — N being the matrix size. However, 
theories for other values of 7] were also considered Q , and the answer for the effective action 
in the leading order in N was obtained in for general 77. 

In the present paper, we derive the constraint equations for the NBI matrix model, 
and show their coincidence with the constraint equations for the Kontsevich matrix model, 
thereby proving the equivalence of the two models for nonzero rj. 



2 The matrix model in the large N limit 

We start with the following matrix integral: 

Z = JdX e -^[* A +^ 1 +(^+ 1 ) lo s*], (21) 

As shown in 0, it is related to the bosonic part of ( |1 . 1| ) by the following change in integration 
variables: Y = ^A _1 , and by 

A = -^[^,^] 2 - (2-2) 

The matrix integral (|2.1| ) belongs to a class of generalized Kontsevich models (GKM) ||. 
Such models with negative powers of the matrix X have been previously discussed in the 
context of c = 1 bosonic string theory 0. In 0, the r-function approach to such models was 
developed. There, the parameter 7] plays the role of the zeroth time in the corresponding 
integrable hierarchy. Moreover, at the conformal point 77 = 0, this model was shown |7j to 
have the same Schwinger-Dyson equations as the U(N) model solved in [|l(], 11]. 

For the models of this type, the large iV solutions are known explicitly only in some special 



cases. The models with cubic potential for X [12] and the combination of the logarithmic 



and quadratic potentials |Qj were solved by a method based on Schwinger-Dyson equations, 



developed first for the unitary matrix models with external field fT0| , |TTf . The same technique, 
being applied to the integral ( j2.1|) , also allows one to find its large iV asymptotic expansions 
in the closed form for arbitrary 77 . 
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The Schwinger-Dyson equations for ( |2.1| ) follow from the identity 
1 d d r .„ d 



N 3 dA ik dA H 



J dX JL e -iVtr[XA + X- + (2, + l)logX] = Q ^ 



Written in terms of the eigenvalues, these N equations read (no summation over i is 
implied) 



I d 2 I I ( d d \ I d 



Z(X) = 0. (2.4) 



For T) — 0, these formulas coincide with the corresponding formulas for the U(N) model [ID 
0- 

It is convenient to set 

This quantity plays an important role in evaluating the large N limit. 
We also introduce the eigenvalue density of the matrix A: 

p(x) = ^J28(x-X i ). (2.6) 

i 

The density obeys the normalization condition 

dxp(x) = l, (2.7) 



and in the large iV limit it becomes a smooth function. 

Simple power counting shows that the derivative of W(Xi) in the first term on the left- 
hand side of Eq. ( j2.4|) is suppressed by the factor 1/N and can be omitted at iV = oo. The 
remaining terms are as follows: 



- xW 2 {x) - I dypiy) y — ^ + (2 V - l)W(x) + 1 = 0, (2.8) 

J y — x 

where Aj is replaced by x. Equation ( |2.8[ ) can be simplified by the substitution 

W(x) =xW(x) -rj. (2.9) 
After some transformations, using the normalization condition ( p.7| ), we obtain 

W\x) + x [ dyp(y) W ^ - W ^ =x + rf. (2.10) 
J y — x 



The nonlinear integral equation (|2.10|) can be solved with the help of the ansatz 



where f(x) is an unknown function to be determined by substituting fl2.11p into Eq. ( |2.10| ). 
The asymptotic behaviors of W(x) and f(x) as x — > oo follow from eq. (|2.10|) : W(x) ~ 
tJx — 1/2, and the analytic solution with minimal set of singularities is 



f[x) = yax + b. 



(2.12) 



The parameters a and b are unambiguously determined from Eq. ( |2.10| ) . We find that b = rj 2 
and a is implicitly defined by 



dy 



p{y) 
f(y) 



or, in terms of the eigenvalues, 



1 + 



2N 



E 



i 



/ \/a\j + 7] 2 

Then, the answer for the integral in the large N limit reads 



(2.13) 



(2.14) 



logZ 



N 2 



2 1 \ Ar] 2 i] 2 

4/ Ja a 



v a\i + rf — rj 

a\i + rj 1 + 7] log I Xi 
« 1 \ -laXi + rf+r] 



~2 E lo s ( v a ^ + ^ 2 + v a ^ + ^ 2 ) • 



(2.15) 



One can verify directly that 4- log Z = and 4 log Z = W(Xi), as far as Eq. ( p.!4| ) holds. 



3 The Kontsevich phase 



We are interested in the asymptotic expansion of the model fl2.1|) for large A. Then, the 
expansion parameters are traces of negative powers of the external matrix A. Conventionally, 
this regime is called the Kontsevich phase of the solution. 

Here an important note is in order. In , we did not discuss which branch of the root- 
positive or negative — should be chosen in ( |2.14j) , since both choices led to the same answer 
for the integral in the large limit ( 2.15|) . However, in what follows, we must fix this sign. 

The Kontsevich phase is the strong coupling regime where the expansion in negative 
powers of Aj is to be performed. Then we see that the dependence is only on Aj n 1 ^ 2 , 

71 = 0,1,.... 

As the first step, we perform the phase analysis for the toy case where all Aj's coincide, 
so ( |2.14j) becomes 

n a r , 2 = 4(aA + rj 2 ), a > 0, A > 0. (3.1) 
(1 - v/a) 2 

Then, obviously, the sign of the square term in ( |2.14j ) is negative for a > 1 and positive for 
< a < 1. 
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Algbebraically, there always (except if 77 = 0) exist two solutions to (|3.1|): one with 
< a < 1 and another with a > 1. For the Kontsevich phase to be possible, we demand 
that the expansion in terms of the so-called times 



1 



2k - 1 

should make sense. So, we assume that 



trA -fc+i/2 = 1,2, 



< 1. 



Then, from (|3.1| ) and (|3.3|) , we obtain the following phase diagram: 



(3.2) 



(3.3) 



1 

2V2 



no Kontsevich^ 
phase / 3 Kontsevich 
phase for 
y/a > 1 

two 
Kontsevich 

''phases for y/a > 1 

and \fa < 1 



1 

2V2 



4 Constraint equations in the Kontsevich phase 



Now we write ( |2.4| ) in terms of the relevant times ( |3.2j ). Here, to obtain rigorous results, the 
normalizing factor is necessary. From the theory of the generalized Kontsevich model ||, 
the proper expression, which has no explicit dependence on the matrix size N, reads 

/ DXe AX+v W 
Z({T " }) = e tr A^o+v(x ) de t-V2 - -S_ V (x )) ' 

Here Xq is the stationary point, A + V'(Xq) = 0, and the determinant in the normalizing 
factor comes from the quasi-classical integration. 
Let us choose the new variables 

A* = zt - (rj + 1/2) 2 . (4.2) 

Then the normalizing factor reads 

exp j-iV^pz, - 2r ] \og(z l + 77 + 1/2)] - ± £ log(z, + zA . (4.3) 



Note that if we perform the standard Itzykson-Zuber integration, expression (JO) for 
integral (|2.1| ) becomes 



Z{{r n }) 



det 

l<i<jV 
0<i<iV-l 



2rjN-l 



2riN-l 



rii e- 2N ^{z i + r ] + l/2)^ N {2z l )- 1 / 2 A{z) 



(4.4) 



where K v (x) are Macdonald functions, K v (x) = J^ds e~ 2xcoshs+lJS , and ( |4.4j ) does not 
resemble too much the corresponding expression for the Kontsevich matrix model where 
these Airy functions stand instead of K u (x). To find the large N asymptotic behavior of 
we use the constraint equation method. 
In terms of ^-variables, fl2.4|) becomes 



N 2 



2z; dzj \ 2z; dzj 



1 zj-(r? + l/2) 2 
A?" 2 2 2 - z 2 



1 9 



1 9 



2zjdzj 2zidzi 



277-I 1 9 

iV 2^9^ 



+ 1 



(4.5) 



When pushing the normalizing factor ( (4.3| ) through derivatives w.r.t. ^-variables, we replace 

1 ^1 



9i -> 9j - 2iV + 2r/iV 



*i + T] + 1/2 



E 



Zi -\- Zj 



Let us introduce the new times 

1 



2^T?^ + ^°^Ti7^' n = ' 1 --- 



(4.6) 



(4.7) 



which differs slightly from the conventional ones defines above by ( |3.2| ). As is easily checked, 
they are related by a lower triangular transformation, so they are equivalent from the view 
point of phase transitions and critical behavior. 

Then, the constraint equations for Z({t n }) are obtained after some tedious algebra which 
we omit here. Collecting all coefficients to the term = -^-t s+ 2, we obtain 

z- azi 



L s Z({t n }) = 0, 



where 



16iV 2 



+ tl 



(?7 + l/2) ; 
4iV 2 



s > -1, 



16N 2 



(4.8) 



2iV 2 



1 oo 

— £ [-(r/ + l/2) 2 (2p+l)t p+1 + (2p-l)(l-5 Pi0 )t ? 



p=0 
p+s>0 



d 



Of 



p+s+l 



N \dt s+2 

1 'f 9 



at, 



+i 



5 



| (rj + 1/2) 2 ^ f 7 ^ 



4AT2 



k=l 



dtk dt 



s+l-k 



(4.9) 



5 



In the convenient conformal field theory notation, 



becomes 



L s = K+i - (V + 1/2) V. + 



7] + 1/2 



«2s+3- 



where 



«2n+l 

n>0 

[ a -2k-l, «2<j+l 







= (271+1)^1, 

-(2k + l)5 k , q , 



and the operators V s in the free- field representation are quadratic in a*, 

^9 = X!^a+6+i:a2a+ia26+l : + l/45 s ,o , 
tt.6 



(4.10) 

(4.11) 

(4.12) 
(4.13) 

(4.14) 



where the normal ordering : : means that all a a with positive indices are on the right of all 
ab with negative indices. 

For s,t > — 1, L s satisfy the algebra 



[L s , L t ) = 4(s - t)(L s+t+1 -(rj+ 1/2) 2 L 



S+t), 



(4.15) 



from which the Virasoro algebra can be obtained by the lower-triangle replacement 

oo 

L s = L s + + l/2)' 2k L s+k . (4.16) 

k=l 

Performing replacement (|4.16|) and rescaling V s , we obtain the Virasoro algebra in terms of 
the generators 

OO 1 

L s = -( V + 1/2) V s + A V N g {7] + 1/2)2k . 3 d s+k . (4.17) 

Amazingly, if we manage to remove the derivative terms in ( |4.17| ), then the constraints we 
obtain will be just constraints of the Kontsevich matrix model [|13| , [L4]. To remove these 
terms, we shift all of the higher times, leaving the times t\ and t% unshifted, 



tk = th — 



4t]N 



( v + l/2) 2k + 1 2(2k + l)' 
Explicitly, in terms of the new times, 



k > 3, 



U 2 = t 



1,2- 



00 d AnN 

L_ 1 = P l + 2E(am + l)Wx K - 5n L ? ft 



(4.18) 



(4.19) 



m=l 

* d d ~ N ~ d 1 AnN d 

■ = E + 2 ,S (2m - "'-c; + ^ " STWi^' (4 - 20) 



This is nothing but the Virasoro algebra that appears in the Kontsevich matrix model. 
Therefore, we have proven the equivalence between the model ( |2.1| ) and the Kontsevich 
matrix model 
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5 A large N limits comparison 

Let us explicitly compare the model (|2.1|) after the time changing fl4.18|) and the Kontsevich 
model with the partition function 



Z(M) 



J dX exp - tr 7 {MX 2 /2 - iX 3 /6} 



JdXexp-tr 1 MX 2 /2 
Let us consider the constraint equation ( |2.14j ). We have 



M— diag{mi, . . . ,m N }. 



(5.1) 



1 



2n- 1 



trM 



-2n+l 



1 



2n - 1 



trZ 



-2n+l 



+ < 



N 



17+1/2' 

0, 



1 



4r)N 



i.e., 



trM -2n+l = tlZ -2n+l 



1 



4r]N 



2(2n-l) (rj+1/2) 2 "- 1 ' 
1 4r]N 



n—1, 
n = 2, 
n > 3, 



+ 



5 n o + 2N8 n s- 



(5.2) 



(5.3) 



2( V + l/2) 2n - 1 2 (77 + 1/2) 3 ™ 

Then we obtain, in terms of these shifted variables, (± depends on the branch of the square 
root), 



1 ± 



1 

2N 



E 



1 



az] - 0(77 + 1/2) 2 + 

00 



I 00 (2n - IV I 

..2n+l 

.(( 7? + l/2) 2 -^/o) n 



2iVv^Y^o (2n)!! 



i±o^ E 



2A^Va 



^m] + 7] 2 /a -{r] + l/2f 



+2N^-2N 
Assuming the minus sign and denoting 

(77 + 1/2) 2 - rf/a = 2s 
we obtain 



7]N 



+ 



r] 3 N 



r] + 1/2 (r/ + l/2) 3 a 
2rjN 



N 

E 



1 



1 \ mj — 2s 



(V + 1/2) 3 



-7s, 



= 7, 



(5.4) 

(5.5) 
(5.6) 



i.e., the constraint equation of the Kontsevich model itself |T3|, TA . 

Also, let us recall the answer in the large N limit for the Kontsevich model |L3]. If 
we reconstruct the dependence on the coupling constant 7, then it reads (in the original 
notation) 

F = 1 2S j + lj:[ m 5 -( m 2 -2sf 2 -3s^^) 



1 Jm 2 -2s + Jm 2 - 2s 

~o E lo s — ! 

2 rrii + raj 



(5.7) 
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where s is determined from the constraint equation 

Let us compare the answer (|5.7| ) with formula ( 2.15Q while accounting for the normalizing 
condition Then, in variables Zi, assuming the minus sign in the constraint (|2.14|) , we 

have 



F = N 2 rf log a + 
+ ^£2 



2^2 



AN r] N 2 r] 



2s + Nrj J2 log 



,(z, + r/ + l/2) 2 



2^ l0g 

*,3 



'zf - 2s + Jzj - 2s 



Zi -\- Zj 



z 2 - 2s 



v/y/a \ 



z 2 - 2s 



(5.i 



After a tedious algebra, taking into account (|5.3| ), we obtain 



F = F + 2N 2 rf log 



V 



V + l/2 



+ N 



7 2 

-V +V 



(5.9) 



where F is given by ( |5.7|) with the substitution ( |5.5|) . The difference between F and F 
is just irrelevant constant terms. This again proves that the two matrix models under 
consideration coincide. 



6 Higher genus expressions 

In this section, we set the Kontsevich coupling constant 7 = 1. The higher genus contribu- 
tions in the Kontsevich model are expressed in terms of the so-called moments Ik, 

1 ^ 1 



Then, for g > 1, we have [14 



1 jh jh r'3 9 -2 

t 9~ 2^ \ T 2 T 3 T 3g-2/g^ r ^2f„-lUVl. LI LI L _ I ' ^ 



39-2: 



i.e., F 9 is a finite sum of monomials in Jfe/(1 — J 1 )( 2A -'+ 1 )/ 3 with coefficients being the inter- 
section indices on the moduli space [pL5f| . For g = 1, Fi = ^ log 

Now we rewrite the expression ( |6.2j ) via the moments of the model ( |2.1| ). Let us introduce 
the new moments Jk 



1 w 



S 



We are interested only in transformation law for with k > since the only dependence on 
the moment Jq is via the constraint equation ( |5.6|) . Then, Eq. ( |5.3|) implies 



i.e., the expression (p.2|) becomes 



(Ii-l)^a 3 / 2 (Ji+W)> 
J fc ^a fc+1 / 2 (J fc + 2/r / 2fc ), 



-, 3g-2 ( J _|_ 2 \ifc 

F NBI _ / I 2 1 3 . . . *3 9 -2v f TT ^ k I ^ > 1 

£ fc=2 (fc-l)**=3 S -3 ( J l + ^) ^ *==2 

(6.4) 

and 

F 1 ^ = llog[a 3 / 2 (J 1 + 2/r / 2 )]. (6.5) 
Therefore, expression ( [2.15 ) for genus zero, taking into account the normalizing factor 



( f4.3|) and the expressions ( |6.5|) , (|6.4|) completely determine the partition function of the 
model ( |2.1| ) for all genera. These expansions are, however, ill-defined for 77 — > 0, which 
corresponds to the U(N) model, and for 77 ~ 1/N, which corresponds to the model defined 
in 0. 



7 Remarks 



1. The last observation above is related to the initial constraints ( |4.8| ). Note that we can 
consider a "minimal reduction," where all times tk but the time t\ are equal to zero. Then 
this partition function is entirely determined from the action of L_i on Z, 



where 



i.e., 



v + 1/2 + ztl at! 



- N 2 (r] + l/2ytj + 1/4 = 



t 2 =t 3 = 



N 



N 



-ii 



l*2=*3= 



^E^r 1/2 + (v +1/2)- 1 , 
iv i=i 



V + 1/2 



iV 2 ((77+ 1/2)^- 277) + 



t2=*3=.»=0 

Remark that T becomes a polynomial in t\ for 



(27/iV) 2 - 1/4 
(77 + 1/2)*! + 277 



(7.1) 



'4JV' 



and this is exactly the point corresponding to the IIB superstring model ||. 

On the other hand, from expressions fl6.4|) , fl6.5|) it is clear that in this case all Jfc 





for k > 0, so there is no ti-dependent contributions from g > 1 and only log a comes 



from the genus one term. One can compare with expression for F\ coming from (|7.1| ) and 
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taking into account the constraint equation ( |2.14j ), one finds an exact coincidence of the two 
expressions, i.e., 

F 1 | t2=i3= ... =0 = -1/8 log(2 + l/N tr A~ 1/2 ) = log a. 

16 

2. The string susceptibility in the large N limit can be obtained by differentiating twice 
w.r.t. the string coupling constant rj. One should check that after the first differentiation of 
fl2.15|) , the stationary condition still holds, so the total derivative coincides with the partial 
one and, as a result, we have 

d 2 logZ dHogZ „ 



dr) 2 drf 



2iV 2 (loga + 3). (7.2) 
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